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We present a model for a class of random binary lattices by introducing a one-dimensional system 
where impurities are placed in one sublattice while host atoms lie on the other sublattice. The 
source of disorder is the stochastic fluctuation of the impurity energy from site to site. We study 
the optical absorption spectra and the peculiarities of the motional narrowing effect at the band 
edges for perturbative and nonperturbative degrees of disorder. Analytical results agree well with 
numerical simulations. 



I. INTRODUCTION 

Several years ago, Knapp raised the concept of mo- 
tional or exchange narrowing in one-dimensional (ID) 
disordered systems [Q] and fruitfully explained many op- 
tical phenomena in low-dimensional disordered systems 
like J aggregates and conjugated polymers (for a review, 
see Refs. and references therein). Roughly speak- 
ing, this author demonstrated that the disorder seen by 
the quasi-particles of the excited system is reduced as 
compared to the seeding degree of disorder, namely the 
width of the probability distribution of the on-site en- 
ergy, as soon as the states of the individual molecules are 
collectivized due to the intermolecular interactions. The 
suppression factor depends on whether or not the disor- 
der is small enough to be regarded as a perturbation. In 
the perturbative case, this factor is determined by the 
square root of the number of sites in the entire system, 
while in the nonperturbative case this number should be 
substituted by the so-called number of coherently bound 
molecules, due to the occurrence of the Anderson local- 
ization whenever the on-site energy is an uncorrelated 
random variable. 

Numerically simulated absorption spectrum of poly- 
silane with an uncorrelated Gaussian distribution of 
nearest-neighbor couplings are similar to those for an un- 
correlated diagonal disorder Q]. In contrast, simulations 
of off-diagonal disorder given by Gaussian randomness 
in the molecular positions |^ found that the behavior 
of the optical observables does not fit the standard mo- 
tional narrowing effect [Q . Recently, this discrepancy has 
been uncovered by Malyshev and Dominguez- Adame , 
pointing out the appearance of correlations in the inter- 
site couplings even if the fluctuations of the molecular 
positions are uncorrelated. Furthermore, it is nowadays 
well established both theoretically and experimen- 
tally 1^ that extended states may arise in ID random sys- 
tems when disorder is correlated, in contrast to the earlier 
believe that all eigenstates should be localized [|l^. The 
motional narrowing effect and the competition between 
the localization length and the correlation length have 



been considered very recently by Rodriguez et al. pT) . 

In this work, we report further progress along the lines 
in the preceding paragraphs. In particular, we focus our 
attention on a new ID tight-binding model of correlated 
disordered system supporting extended states |p^ , |l3| . 
We built up our model by considering a ID random bi- 
nary lattice with two species, referred to as A and B 
atoms hereafter. We further assume that the site en- 
ergy of A atoms is randomly distributed from site to site 
while that of B atoms is the same over the other sublat- 
tice. We have demonstrated analytically the occurrence 
of extended states in the vicinity of the site energy of B 
atoms in spite of the fact that the system is purely ID and 
random p3| . Here we study the peculiarities of the mo- 
tional narrowing effect and their manifestation through 
the optical properties of this system. 



II. MODEL 

Let us consider a tight-binding Hamiltonian with 
nearest- neighbor interactions 
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where the state vector |n) represents an excitation at 
site n. In the present ID binary system, A (B) atoms are 
placed at odd (even) positions of the otherwise regular 
lattice, whose corresponding site energies are e2n-i (£2ri) 
with n = 1,2, ... N, N being the number of unit cells of 
the system. According to our model, site energies at even 
positions are the same and we can set e2n — without 
loss of generality. The source of disorder arise from the 
stochastic fluctuations of site energy at odd positions. 
We assume that {e2ri-i}^=i is a set of uncorrelated ran- 
dom Gaussian variables with mean value v and variance 
a^. Hereafter a will be referred to as degree of disor- 
der. The joint distribution function is represented by the 
direct product of single Gaussians. Thus 
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where the angular brackets (. . .)av indicate the aver- 
age over the ensemble. Although the system is one- 
dimensional and random, it has been demonstrated an- 
alytically the existence of a delocalized state in infinite 
systems ai E = Q |]l3|. Most important, there exist many 
{'-^ ^/N) states close to the resonant energy that remain 
extended in finite systems, in the sense that their local- 
ization length is larger than the system size. 



III. PERTURBATIVE MOTIONAL NARROWING 
EFFECT 

We now calculate the motional narrowing effect in the 
perturbative limit by considering small degree of disor- 
der a. To this end, we start by writing down the eigen- 
states and eigenenergies of the unperturbed Hamiltonian 
(cr = 0), corresponding to a diatomic lattice with on- 
site energies and v in each unit cell. We assume rigid 
boundary conditions and f > 0. The last restriction 
is only for the sake of clarity in the exposition of re- 
sults, as they are trivially generalized to the case v < 0. 
There are two allowed bands, the first one ranging from 
v/2— (v/2Y + 4J^ up to and the second one ranging 
from V up to v/2 + ^J{v/2)'^ + A,P for v > Q. The nor- 
malized eigenstates of the unperturbed Hamiltonian can 
be written as 
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where K = kTT/{2N + 1) with fc = 1, 2, . . . , and 
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Here the labels ± refer to the upper and the lower 
(— ) band, respectively. The dispersion relation within 
each band is given by 
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The perturbative part of the Hamiltonian, Tip, is di- 
agonal in the site representation and comes from the 
stochastic fluctuations of the on-site energy at odd sites 
around the mean value v. Therefore Tip = y2n=i {i2n-i — 
v)\2n—l){2n—l\. In the if -representation W), the matrix 



elements of Tip are expressed through linear combinations 
of Gaussian variables e2n-i ~ v with zero mean 
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where £ and i' indicate the band. Consequently, they 
also have a joint Gaussian distribution with zero mean. 
In what follows, we will deal with the fluctuations of the 
square of the matrix elements given in (^) since they 
determine the linear optical response of the system. 

Concerning optical transitions, since the states with 
fc = 1 (one per band) carry almost the entire oscilla- 
tor strength of the system, in the perturbative limit the 
optical absorption spectrum is dominated by two Gaus- 
sian peaks located at the bottom (top) of the lower (up- 
per) band. It is easy to demonstrate that the oscilla- 
tor strength is proportional to =F (3^, where and 
(3^ are given by (^, so that the lower energy transition 
becomes more intense. Therefore, we restrict ourselves 
to this transition, whose standard deviation is given by 
- fill with Bl^, = {{K -\np\K',-)^),,. After per- 
forming the summations (o) one can obtain Bkk' ■ Close 
to the zone center, which is responsible of the linear op- 
tical response of the system, the calculation yields 
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where in (Q) we have replaced Ej^ and E'^, by their values 
at the zone center since K ^ 0. Therefore, the standard 
deviation of the main absorption line is given by 
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As it can be seen from (||), the standard deviation of 
the Gaussian peak scales as N^^^"^, showing the so-called 
motional narrowing effect Notice that the reduction 
factor carries information on the system parameter v/ J. 



IV. NONPERTURBATIVE DEGREE OF 
DISORDER 

The perturbative approach we have carried out in the 
preceding section holds for small degree of disorder. On 
increasing the degree of disorder, mixing of the states k — 
1 with the other states strongly affect the system optical 
response. Hence, one should compare the energy differ- 
ence between the two lowest states, SE12 = — E]^J, 
with the typical fluctuation of {Ki, — \Tlp\K2, ~) repre- 
sented by B12, where Km = m7r/(2A -f 1). The pertur- 
bative approach is valid provided SE12 > B12 and fails 
otherwise. The energy difference can be readily deter- 
mined from (|^) 
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For higher values of the degree of disorder {SE12 < B12) 
the perturbative approach fails; not all sites contribute 
to the optical spectrum since the Anderson localization 
length becomes smaller than the system size. Thus, 
N should be replaced by N* < iV in the above equa- 
tions 01 , where N* is often referred to as the number of 
coherently bound molecules. Malyshev proposed a self- 
consistent estimation of N* by applying the perturbative 
criterion mentioned above to a typical localization seg- 
ment of length N* [Q. This approach yields excellent 
results for uncorrelated and pairwise correlated diagonal 



disorder 1^151. Thus, we use the condition B 
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where the asterisk means the substitution of N by N* 
in (ItI) and (p|). In doing so we obtain 
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Finally, replacing in (g) N by N* given by (|10|) we get 
for the standard deviation of the peaks 
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where 



C = 



V3 



4/3 



2(12\/27r2)i/3 



^(w/ J)2 + 16 



(11a) 



(lib) 



This expression holds for w > as well as for z; < 0. 
The power-like dependence a^l'^ also appears in dealing 
with uncorrelated and pairwise correlated diagonal dis- 
order but with different C (see, e.g., Ref |l^) and can 
be obtained from a more rigorous framework using the 
coherent potential approximation |16|. 



V. NUMERICAL SIMULATIONS 

To check the accuracy of the above analytical ap- 
proach, we have obtained numerically the absorption line 
shape, 1(E), according to Ref. We fix the value J = 1 
and focus our attention on the degree of disorder cr, rang- 
ing from 0.1 up to 0.9. We have diagonalized the Hamil- 
tonian (^ for chains of iV = 125 (250 sites) with rigid 
boundary conditions. The number of randomly gener- 
ated systems is 10 000 for each value of tr. We show in 
Fig. |l| several examples of the optical absorption spectra 
for various values of the parameters v and cr; only the 
lower energy side of the spectra is shown. For small val- 
ues of a the position of the main line agrees very well with 
the value w/2 - ^(u/2)2 obtained from (||). This 

main line broadens and becomes redshifted on increasing 
the degree of disorder, as expected. 
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FIG. 1. Absorption spectra for one-dimensional random bi- 
nary lattices with J = 1, N = 125 (250 sites) and v = — 1 
(solid lines) and v = +1 (dashed lines). The narrow and broad 
peaks correspond to cr = 0.1 and a = 0.9, respectively. Notice 
the scale factor in the later case. Each spectrum comprises 
the results of 10 000 realizations of the disordered system. 

The value of al was obtained by nonlinear Gaussian 
fitting of the low energy side of the spectra and the results 
are presented in Fig. y. It is found that this magnitude 
scales as a^^^ for the var ious values of v studied, as pre- 
dicted by our estimates (11a). The slopes of the st raigh t 
lines are slightly larger than those obtained from ( |lll:| ), 
as shown in Table |; nevertheless, the coincidence should 
be admitted as being highly surprising in view of the 
simple assumptions we made. 
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FIG. 2. Standard deviation ai as a function of o"*^^ for 
different values of v, indicated at the right of each line, with 
J — 1 and — 125 (250 sites). Solid lines represent the 
least square fits and each point comprises the results of 10 000 
realizations of the disordered system. 
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VI. CONCLUSIONS 

In this paper we have considered the motional narrow- 
ing effect in a ID random binary lattice where disorder 
lies in one of the two sublattices. This system could be 
regarded as a simple model of semiconductor binary alloy 
— like ternary III-V compounds — . In these alloys (say 
A^Gai_2;As), the cation sublattice is occupied by the 
same atoms (say As) while anions (say Al and Ga) are 
randomly distributed over the other sublattice. Starting 
from a pcrturbative approach, we found that the width 
of the absorption line scales as ~ N~^^'^ for small degree 
of disorder, showing the motional narrowing effect. For 
larger degrees of disorder, we determined self-consistently 
the spatial extend of the excitation wave function, ac- 
cording to the prescription raised by Malyshev . The 
spatial extend of the excitation wave function is smaller 
than the system size and the width of the absorption line 
scales as ~ a'^^^. Our estimates agree well with numeri- 
cally simulated optical absorption spectra. 



[13] F. Dommguez-Adame, I. Gomez, A. Avakyan, D. Se- 

drakyan, and A. Sedrakyan, phys. stat. sol. (b) (accepted 

for publication). 
[14] V. A. Malyshev, Opt. Spectrosc. 71, 505 (1991); J. Lu- 

min. 55, 225 (1993). 
[15] F. Dommguez-Adame and V. A. Malyshev, J. Lumin. 

83-84, 61 (1999). 
[16] A. Boukahil and D. L. Huber, J. Lumin. 45, 13 (1990). 



v/J 


-2 


-1 





1 


2 


10 


Analytical 


0.42 


0.34 


0.25 


0.17 


0.12 


0.010 


Numerical 


0.57 


0.51 


0.40 


0.29 


0.19 


0.014 



TA BLE 1. Slope of al-a'^^^ plots obtained analytically 
from (lib) and numerically for different values of the param- 
eter V. 
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